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To understand the electronic shell- and supershell-structure in large metal clusters we have 
performed self-consistent calculations in the homogeneous, spherical jellium model for a variety of 
different materials. A scaling analysis of the results reveals a surprisingly simple dependence of 
the supershells on the jellium density. It is shown how this can be understood in the framework 
of a periodic-orbit-expansion by analytically extending the well-known semiclassical treatment of a 
spherical cavity to more realistic potentials. 
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The structure observed in the mass spectra of large, 
warmed metal clusters can be attributed to the 

properties of itinerant electrons moving in a finite volume 
[§,[3|. The most prominent finite-size effect is the occur- 
rence of pronounced oscillations in the density of states 
JjJ giving rise to an oscillating part E of the total energy, 
which is superimposed on the smooth Thomas-Fermi en- 
ergy E. With increasing cluster radius, one finds regular 
oscillations (shells) whose amplitude is modulated (su- 
pershells) . 

Two different theoretical approaches have been used 
to describe the electronic shell- and supershell-structure 
in large metal clusters. One is the self-consistent jellium 
model ||. In its simplest form a cluster is described by 
a homogeneous sphere of given charge density, dressed 
with TV valence electrons. Treating the electrons self- 
consistently using density functional theory, the only in- 
put parameter for such a calculation is the Wigner-Seitz 
radius r s . Although describing the electronic structure of 
alkali clusters quite well, this model provides little physi- 
cal insight into the mechanisms underlying the shell- and 
supershell oscillations. It is here that the second ap- 
proach comes in. Given an effective one-particle potential 
one can find a semiclassical expansion of the oscillating 
part of the density of states in terms of classical periodic 
orbits j7|,[|. Introducing a suitable damping factor one 
finds for the spherical cavity, that the oscillations in the 
density of states are essentially determined by the contri- 
butions of triangular and square orbits. The supershells 
can thus be understood as a beating pattern originating 
from the contributions of these orbits. 

The semiclassical approach draws its power from the 
fact that the periodic-orbit-expansion is known analyti- 
cally for the model case of the spherical cavity. However, 
potentials with hard walls are only a crude approxima- 
tion to realistic cluster potentials, which have a soft sur- 
face. Such potentials can also be treated using semiclassi- 
cal techniques Jl(],[ll| ■ In particular ultra-soft potentials 
have received much attention [ ^2p^ |. Unfortunately in 
these cases the action integrals entering the semiclassical 
formalism have to be evaluated numerically. 
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In order to find out what determines the shell- and 
supershell-structure in metal clusters, we have carried 
out the following program. First we performed a series of 
calculations in the homogeneous, spherical jellium model 
for a range of different electron densities. A scaling anal- 
ysis of the results suggested that changing the electron 
density merely introduces a phase shift in the supershell 
structure. It is shown how this can be understood semi- 
classically in terms of a leptodermous expansion, where 
the action integrals for the potential under consideration 
are expanded around a cavity. Thereby we obtain an 
analytical expression for the shift of the supershells. Fi- 
nally it is shown that the leptodermous expansion works 
for realistic cluster potentials by comparing the shifts of 
supershells extracted from self-consistent calculations to 
those given by the semiclassical formulae. 

As starting point of our analysis, we have performed 
extensive calculations in the homogeneous, spherical jel- 
lium model. We use the local density approximation 
in the parameterization given in Ref. \L^ . Electron 
densities range from r s = 2.07 a for aluminum to 
r s = 5.63 do, corresponding to bulk cesium. Cluster sizes 
were chosen from N — 100 to 6000 valence-electrons, thus 
including the first two nodes of the supershell oscillation 
for all densities considered. Typical results for E(N) are 
shown in fig. [|. It can be seen that the supershells are 
shifted towards larger N as the Wigner-Seitz radius de- 
creases, while the positions of the shell-minima are fairly 
independent of r s . To quantify the supershells we de- 
termine the envelope of E(N) by low-pass filtering its 
absolute value. The position of the supershell-nodes are 
given by the minima in the envelope. Filtering out the 
shell structure of course introduces an uncertainty of the 
order of the distance between adjacent shell- minima. The 
results of our jellium calculations are listed in table Q. 

To make a quantitative comparison of the different jel- 
lium calculations, we describe the problem in terms of 
dimensionless quantities. In order to do so, the relevant 
scales of the problem have to be identified. Obviously 
one such scale is the Wigner-Seitz radius. In fact, we 
find that the amplitude of the oscillations E(N) are 



1 



11U 




FIG. 1. Oscillating part of the total energy extracted 
from self-consistent calculations in the homogeneous, spheri- 
cal jellium model. The positions of the super-nodes are indi- 
cated. 

proportional to The existence of a surface intro- 

duces an additional scale: the width a of the surface re- 
gion. As has been shown in |l5|] , a is fairly independent of 
r s . Assuming the shift of the supershells to be a surface- 
effect, we identify a/r s oc l/r s as the relevant scaling 
parameter. By plotting the positions of the super-nodes 
as a function of 1 /r s (cf . fig. |2|) we indeed find a simple re- 
lation: The super-nodes are linearly shifted as a function 
of l/r s . In particular, the first and second super-node 
are shifted in parallel. Describing E(N) semiclassically 
as a simple beating pattern Q, it is therefore tempting 
to conclude that the shift of the supershells is caused by 
a phase-shift in the contributions of the periodic orbits. 



material 


r s in a 


1st super-node 


2nd super-node 


Cs 


5.63 


8.39 


14.59 


Rb 


5.20 


8.43 


14.67 


K 


4.86 


8.47 


14.75 


Na 


3.93 


8.95 


15.11 


Li 


3.26 


9.15 


15.45 


Tl 


2.48 


9.79 


16.01 


In 


2.41 


9.85 


16.11 


Ga 


2.19 


10.13 


16.37 


Al 


2.07 


10.20 


16.49 



TABLE I. Position of super-nodes (given as N 1/3 ) for dif- 
ferent jellium-densities. 



To check this conjecture, we derive an explicit periodic- 
orbit-expansion for the oscillating part E(N) of the total 
energy. Our approach is based on the fact that E can, to 
first order, be extracted from the spectrum of smooth po- 
tentials that fit the self-consistent results |l6[ ]. We start 
from the observation jl7) that in the semiclassical ap- 
proximation the density of states naturally separates into 
two contributions: the smooth Thomas-Fermi term p and 
an oscillating contribution p. The latter term describes 
the quantum corrections to Thomas- Fermi theory and is 
given by an expansion over all classical periodic orbits 
which exist in the potential under consideration. For a 
spherically symmetric potential with exactly two classical 
turning points the periodic orbits can be uniquely labeled 
by two positive integers: the number of times A it turns 
around the origin, and the number v of vertices it has. 
Denoting the classical action along such a periodic orbit 
by S{\ %v \i one finds an expression of the form 0: 



■ (S(\, u )/h - ip(\, v )) dE, (1) 



p(E)dE=J2 A M 



where <P(\, V ) is the so called Maslov-phase. Unfortunately 
the expansion (|l|) converges quite slowly. This is obvious, 
since it is supposed to approximate the density of states, 
which is a sum of 5-functions. Therefore one usually in- 
troduces some damping as to broaden the eigenstates and 
make the expansion (|l|) converge more rapidly. But we 
are actually not interested in the density of states, rather 
in E(N). In the limit of large N (which corresponds to 
the semiclassical limit) the oscillating part of the total 
energy is given by 



E(N) 



E F (N) 



dE / dE'p{N;E'). 



(2) 
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2nd super-node 
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FIG. 2. Position of the first and second super-node as a 
function of l/r s . The error bars are due to the uncertainty in 
locating the super-nodes as minima of the envelope of E(N). 
The solid line gives a linear fit to the data, the parameters of 
which are given in the plot. 
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Here E F (N) is the Fermi energy in Thomas Fermi ap- 
proximation. Integrating twice by parts over p essen- 
tially divides the expansion parameters A^\, v ) in (|l|) by 
the square of the classical action along the orbit, thereby 
reducing the importance of the longer orbits: 



E(N) 



El 
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(S(x,v)/h-<P(\,v)y (3) 



Thus there is no need of introducing a damping factor 
to accelerate convergence. Actually, E(N) is dominated 
by the contributions of the shortest plane periodic orbits, 
namely the triangular and the square orbit. Furthermore 
inspection of (||) shows that variations in the boundary 
conditions, which strongly shift the oscillations in p [ p"8[ , 
hardly influence E, since the changes in the density of 
states are compensated by those in the Fermi energy. 

As an immediate application of the expansion (^|) , we 
can investigate how an increase in density for small clus- 
ters compared to the bulk affects the electronic shells 
and supcrshells. Such a lattice contraction was suggested 
by EXAFS analyses of small clusters jl9| . Changing the 
density clearly will change the Fermi-energy for a given 
cluster. As we can see from equation (||) this will obvi- 
ously change the overall amplitude of E(N). But apart 
from that the oscillations are determined by the classical 
actions S(\ ^ along the periodic orbits. For a spherical 
cavity of radius Rq we find 

S(\^)/h = 2v sin(Tr\/v) k F R ; k F — \/2mE F /h, (4) 

where the product k F Ro depends on the number N of 
electrons inside the cavity, but is independent of the 
electron-density. Hence the electronic shell-structure for 
spherical-cavity-clusters does not depend on any lattice 
contraction, except for an overall change in amplitude. 
This result suggests that the same is true for smooth 
potentials, provided the lattice contraction is not too 
large. We have confirmed this by numerically solving 
the quantum mechanical problem for realistic potentials 
introducing contractions ARq of up to 1/2 r s . Thus we 
can conclude that a possible lattice contraction will not 
noticeably affect the electronic shells and supershells. 

Next we turn to the problem of understanding why 
the super- nodes are phase-shifted as a function of jellium 
density (cf. fig. ||). The most straightforward approach 
would be to solve the integrals, which enter equation (j^) , 
explicitly. Unfortunately, this cannot be done analyti- 
cally. But in the semiclassical limit, which corresponds 
to N ^oo, it is sufficient to know the integrals to leading 
order in 1/N. The basic idea is then to use the spherical 
cavity as a starting point and expand the action for more 
realistic potentials around this case, using the surface- 
width a as small parameter (leptodermous expansion). 
For the classical action we then find 



5/t r< cav ity /t 



h +h—\ +0{N- 1 / 3 ), 



(5) 



where the expansion parameters I\ and I2 are indepen- 
dent of cluster-size; i.e. to leading order surface softness 
introduces a phase-shift in the periodic-orbit-expansion 
(||) , while the period of the oscillations is still determined 
by STffl/h. An additional phase shift arises from the 
difference of the Maslov phases f(\,v) for a soft potential 
or a cavity. Finally an inspection of the amplitudes in 
(Q) shows that, to leading order in 1/N, they do not de- 
pendent on the shape of the potential. We thus find that 
replacing a cavity by a soft potential with small surface- 
parameter a amounts to merely shifting phases in the 
periodic-orbit-expansion of E(N). 

Now the question arises, whether typical cluster po- 
tentials are such that their surface-parameter is small 
enough for the above expansion to be valid. To judge 
this, we have to fit the potentials obtained from our self- 
consistent calculation with some analytical model poten- 
tial. Since the classical action depends only on the poten- 
tial in the classically allowed region, it seems reasonable 
to fit only for E < E F . There the self-consistent poten- 
tial, except for possible Friedel oscillations, can be well 
described by a Woods-Saxon function 



V(r) 



-Vo 



l + exp(l^) 



(6) 



But fitting only for E < E F seems to imply an error 
in calculating the Maslov phases which serve to cap- 
ture the influence of the classically forbidden region. 
From WKB quantization it can be seen that the Maslov 
phases for a separable system are given by the sum of 
the quantum-mechanical scattering-phases at the classi- 
cal turning points. For Woods-Saxon potentials we can 
calculate these analytically. To leading order, they coin- 
cide with the Maslov phases for a square-well potential 
of depth — Vo jig]. Thus, the error in the potential for 
E > E F will not enter the leptodermous expansion. 

Given the potential (||), we have found analytical ex- 
pressions for the parameters I\ and I2 in the expansion 
(H) of the classical action. Introducing the abbreviations 



P = 



E F + Vq /ttA 

and Pl = P cos — 

Vq V v 



the expansion parameters for a given periodic orbit (A, v) 
are 
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Details of the calculation will be published elsewhere [^0| . 

We are now in the position to find a simple estimate for 
the shift of the supershells. Following || we start from 
a drastically simplified version of equation (0): All peri- 
odic orbits except triangular (1,3) and the square (1,4) 
orbits are neglected. Furthermore it is assumed that the 
amplitudes for these orbits are equal. This leaves us with 
an expression of the form 

E = A (cos (A JV* + + cos (/ a N$ +cp 2 )). (7) 

Fitting the self-consistent potentials with (^]), we can 
compare the shift observed in the jellium calculations to 
that determined by the leptodermous expansion using 
the ansatz ([?]). This is shown in fig. @. 
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FIG. 3. Shift of super-nodes (in units of N 1 ^ 3 ) for differ- 
ent jellium clusters. The shift obtained from the periodic orbit 
expansion is compared to that observed in the self-consistent 
calculations. If these shifts agreed perfectly, the points would 
fall on the full line. 

Although we have introduced a number of approxima- 
tions the agreement is remarkable. We can thus conclude 
that the leptodermous expansion applies to typical clus- 
ter potentials. 

To summarize, we have demonstrated how the semi- 
classical description of a spherical cavity can be general- 
ized to describe the electronic supershells of realistic po- 
tentials. Starting from a periodic-orbit-expansion for the 
oscillating part of the total energy, we have shown that 
E(N) is hardly influenced by the Maslov Phases or a lat- 
tice contraction. Introducing a leptodermous expansion 
for the classical action we have established that a soft 
potential gives phase-shifts in the semiclassical expres- 
sion for E(N). We can thus understand the dependence 
of the supershell on the electron density, revealed by a 
scaling analysis of our jellium calculations. 

Moreover the leptodermous expansion can be used to 
analyze how the supershell structure will change if the 
underlying model is changed. Introducing a pseudopo- 



tential, as e.g. in the stabilized jellium model |2l[] . in- 
creases the spill-out of the electrons and leads to a soften- 
ing of the potential at the surface of the cluster. This in- 
duces a shift of the super-nodes towards larger N, which 
can be estimated by the semiclassical technique described 
above. 

Finally, the identification of r s as the typical length 
scale for the supershell-problem suggests a justification 
of the ad-hoc procedure proposed in pj to improve the 
results of jellium calculations for gallium clusters. There 
it was found that the introduction of a non-homogeneous 
jellium background is essential for treating Gajv clusters, 
while alkali clusters are well described by a homogeneous 
jellium. Assuming, that the typical length scale for fea- 
tures in the jellium is the ionic radius r at , while the length 
scale for the electrons is the Wigner-Seitz radius r s , we 
find that the importance of inhomogeneities increases 
with the number of valence electrons Z va i oc (r at /r s ) 3 . 

I am much indebted to O. Gunnarsson for his invalu- 
able advice. Helpful discussions with T. P. Martin and 
M. Brack are gratefully acknowledged. 
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